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Abstract 

In this paper, we focus on the role of power corrections in QCD factorization(QCDF) method 
in charmless two-body nonleptonic B meson decays. We use the ratio of the branching fraction of 
B + — » ir + K*° to that of B° — ► ir~ p + , for which the theoretical uncertainties are greatly reduced, 
to show clearly that the power corrections in charmless B decays are probably large. With other 
similar ratios considered, for example, for the B° — ► K~ p + decay, it is very likely that, among 
various sources of power corrections, annihilation topology plays an indispensable role at least for 
penguin dominated PV channels. We also consider some selective ratios of direct CP asymmetries. 
Among these, we find that, if power corrections other than the chirally enhanced power corrections 
and annihilation topology were negligible, QCDF would predict the direct CP asymmetry of B — > 
7r + 7r _ to be about 3 times larger than that of B — ► ir^K^, with opposite sign. Experimentally 
any significant deviation from this prediction would suggest either new physics or possibly the 
importance of long-distance rescattering effects. 
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I. INTRODUCTION 

With the excellent performance of KEK and SLAC B factories, a great number of charm- 
less B decay channels have been measured with high precision. It is even more exciting 
that, with the expected larger data sample in the near future, CP violations in many de- 
cay channels might be soon reachable. However in most cases, our ignorance of strong 
dynamics stands in the way of either identifying the signal of new physics or extracting the 
fundamental Cabibbo-Kobayashi-Maskawa (CKM) angles from the observables. 

Theorists have tried hard to have a better understanding on the strong dynamics in non- 
leptonic B decays, in which QCD factorization method (QCDF)|l[ is a recent progress. In 
the heavy quark limit, QCDF has shown that the hadronic matrix elements of B decays 
have simpler non-perturbative structures, while the QCD corrections are perturbatively cal- 
culable, at least at one- loop order. However the l/m& power corrections, especially annihi- 
lation contributions and chirally enhanced power corrections 

as 

, are phenomenologically 

important. It is even worse that factorization generally breaks down beyond the leading 
power expansion. As a result, the model-dependent parametrization has to be introduced 
in QCDF to account for these formally power-suppressed contributions. For some decay 
channels, these power corrections could even compete with the leading power contributions 
in some parameter space which makes the reliability of QCDF predictions in doubt. Very re- 
cently, a new effective theory, called soft-collinear effective theory (SCET), has been applied 
to charmless B decays in which the power corrections can be studied in a systematic 
way, though more non-perturbative operators have to be introduced. However since this 
method is still under development and several issues remain to be resolved, we shall not 
discuss SCET further in this paper. 

Thus it is of great interest to investigate the role of power corrections in charmless B 
decays, and to determine if the power corrections were really important, whether they were 
mainly from chirally enhanced corrections and annihilation topology or from other sources 
of power corrections, such as long-distance rescatterings. In this paper, we shall discuss 
these problems in a transparent way based on experimental measurements. 
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II. QCD FACTORIZATION FOR CHARMLESS B DECAYS 



To make the paper self-contained, we will recapitulate the main point of QCDF in the 
following. One may find more details in Refs 

There are three distinct scales: M w 3> m b 3> Aq C d involved in charmless B decays. To 
go beyond naive model estimation, the physics of different scales should be separated from 
each other. It is known that, with the operator product expansion and renormalization group 
equation, the effective Lagrangian can be obtained, in which short- distance effects involving 
large virtual momenta of the loop corrections from the scale My/ down to \i = 0(m&) are 
cleanly integrated into the Wilson coefficients. Then the amplitude for the decay B — > M\M 2 
can be expressed as |z|: 

A(B -> M,M 2 ) = ^=E £ VMpKMiMalQiMlfl) , (1) 

V * i=l q=u,c 

where X q is a CKM factor, Cj(/x) is the Wilson coefficient which is perturbatively calculable 
from first principles. The effective operators Qi,2 and Qz,...^ are tree level and QCD penguin 
operators, respectively. We have neglected electroweak penguin operators <57,...,io because 
their effects are numerically small in most decay channels. {MiMi\Q%{}l)\B) is a hadronic 
matrix element which contains the physics effects from the scale \x = C?(m&) down to Aqcd- 
Since the perturbative and nonperturbative effects related to m& and Aq C d still entangle 
with each other, it is highly nontrivial to estimate the hadronic matrix elements. But 
in the heavy quark limit, QCDF shows that the above hadronic matrix elements can be 
factorized into hard radiative corrections and simpler nonperturbative structures which can 
be parametrized by the form factors and meson light-cone distribution amplitudes (LCDAs). 

But for phenomenological applications, the power corrections in l/m^, especially the 
chirally enhanced corrections and annihilation contributions should be taken into account. 
Chirally enhanced corrections arise from two-body twist-3 LCDAs of the final state mesons. 
Unfortunately the factorization breaks down for hard spectator scattering diagrams because 
the twist-3 LCDA does not approach zero in the endpoint region which introduces a log- 
arithmic divergence. A similar divergence also appears in the annihilation contributions. 
Phenomenologically, Beneke et al. |2j adopted a model parametrization for the endpoint 
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divergence: 

Xa,h = / — = log -^(1 + p A ,ne^ H ) (0 < p A , H < 1) , (2) 

J X I\h 



where {Xh) denotes the divergent terms from annihilation topology (hard spectator 
scattering). The corresponding price to pay is model-dependence uncertainties and also 
large numerical uncertainties from Xa,h- With the above discussions, the decay amplitudes 
can be written as 



A(B - M 1 M 2 ) = ^ £ v p [Y,<(M 1 M 2 \O i \B) f + Y,fBfMjM 2 b 

V * p=n,c \i=l 7 



(3) 



where (M1M2IO1I-B) / is the factorized hadronic matrix element which has the same definition 
as that in the naive factorization approach. For the complete expressions of QCD coefficients 
Oj and annihilation parameters bj, one may refer to Ref. p. 

III. ARE POWER CORRECTIONS LARGE ? 

It is clear that chirally enhanced power corrections should be large, at least in B decays to 
two light pseudoscalar mesons, because although they are suppressed by the factor (taking 
the pion as an illustration) 

2m 2 

r x = 1 \ > ( 4 ) 

* m b (m u + m d ) 

which is formally of the order of Agc^/mf,, numerically this factor is of order one. Experi- 
mentally, the large branching ratios of B — ► tiK decays also strongly support this statement. 
For example, the ratio 

B(B+ -> vr+if ) _ (21.8 ± 1.4) x IP' 6 , v 

B{B+ -> n+K*°) (9.0 ± 1.3) x 10- 6 ' [) 

would be theoretically smaller than 1 without chirally enhanced contributions, say r^a^ 
term. The experimental data in the above equation are from ref. jsj. 

However it is not very clear about the role of annihilation topology in B decays. It was 
once believed to be quite small in most B decay channels, because of the power suppression 
I9JI. The importance of annihilation contributions was first noticed in the jperturbative QCD 

n 



method for charmless B decays [10]. Recent phenomenological analyses p, [ill, |12|, |13J| based 



on QCDF also suggest substantial contributions from annihilation topology. However since 



in QCDF, there are many parameters involved in the global analysis of experimental data, it 
would be interesting to show the importance of annihilation topology in a transparent way. 
In the following we will try to do it with various ratios. We find that annihilation topology 
very likely plays a significant role, at least in penguin dominated B — ► PV decays, where P 
denotes light pseudoscalar meson and V denotes light vector meson. 

Let us first consider the ratio of the decay amplitude of B + — > tt + K*° to that of B° — > 
p + 7r~. If the power corrections were negligible, this ratio would be theoretically very clean 
where the form factors cancel out, furthermore it is almost independent ot the CKM angle 
7 and the strange-quark mass: 



A(B- 



7l + K 



*0\ 



A(B C 



V cb V cs 



VubVud 



fa 
fp 



a c 4 {7iK*)+rfa c 6 {7iK 



n " 

a L 



(6) 



where the penguin contributions to B° — ► p + ir decay and the term proportional to V u bV us 
in the numerator are neglected. This should be a reasonable approximation up to a few 



percent level. In QCDF, \(al(%K*) 



n 

(With the default parameters in ref. [a], it is 0.03) and fx* 
for the CKM matrix elements, there is a useful inequality 



r!s ag(TrK*))/d^\ should be about or less than 0.04 




is very close to unity. While 



V, 



ab 



V, 



cb 



Asin/3W1 + 



cos 2 a 



sin 2 a 



> A sin (3 . 



(7) 



The current measured value of sin 2(3 from BaBar and Belle gives 



sin 2(3 = 0.736 ±0.049 



□ 



(8) 



from which we obtain sin/3 = 0.402 ± 0.033. In principle we could get another solution 
sin (3 ~ 0.9 which however is inconsistent with the direct \V u b\ measurements. So given 
A = 0.224, we will get an interesting lower limit 

Vub 



cb 



> A sin (3 = 0.090 ± 0.007 > 0.078 ( 90% C.L.) . 



(9) 



With this lower limit, there arises a clear discrepancy between theory and experiments: 



0.53 > 



A(B- 



7T 



K 



A(B C 



r{B°) B{B + -> 7i+K*°) 



1/2 



0.77 ±0.09 



(10) 



T(B+) B(B° -> P+7T-) 

where the left hand side is from theoretical estimation Eqs. ©, ©, © and the right hand 
side is obtained from the following experimental measurements: B(B° — > p + 7r~) = 13.9 ±2.7 



13, B(B + -> 7i+K*°) = 9.0 ± 1.3 



and t(B+)/t{B°) = 1.083 ± 0.017 Q- Notice that 
cos 2 a is of the order of a few percent, it seems appropriate to neglect it in Eq. (JJJ). Then 
the left hand side of Eq. (JTUJ) would be reduced further to 0.46 ± 0.04 if we fix (a 4 + r x a 6 ) / a\ 
to be 0.04. In fact, (a 4 + r x ao) / a\ for nK* channels is typically less than 0.04 in most of the 
QCDF parameter space. For example, if using the default parameters in ref. |5j, we will get 
0.35 for this ratio. Since the chirally enhanced corrections for this ratio are not expected to 
be large because ae(7i + K*) vanishes at tree level, this is a strong indication that annihilation 
topology and/or other sources of power corrections might play an important role at least in 
B — > PV decays. 



If we turn to another ratio, the branching fraction of B° — > K + p to that of B c 



p 7T~ 



there is similar disagreement, although this time it is not well established considering the 
large theoretical uncertainties. Assuming negligible annihilation contributions, the form 
factors cancel out again for this ratio. However there are significant dependences on angle 
7, V u b/V c b and especially the strange-quark mass m s in the numerator because of the factor 



A(B° -> K + p 



A(B° -> p-n+) 



f_K 
U 



cb 



V, 



ub 



ai(pK)-rt<a c 6 (pK) 



(11) 



where as usual, the small penguin terms in the denominator has been neglected. If we take 
the default parameters in ref. p, namely 7 = 70°, V u b/V c b = 0.09 and a%{pK) —r^al(pK) = 
0.037 + 0.003z with the strange-quark mass set to be 90 MeV, the theoretical prediction 
for the ratio B(B° K + p~)/B(B° — ► p~7i + ) is 0.38, which is again significantly smaller 



than the experimental measurements 1.01 ± 0.34 |8l [TTf. Of course, when the theoretical 
uncertainties, especially the uncertainty of the strange-quark mass, are taken into account, 
the disagreement is not that impressive. For example, the ratio can rise up to 0.69 when m s 
is lowered to 70 MeV. 

However when combining with the B + — > tt + K*° decay, this result shows that there is 
a tendency that, the penguin-dominant B PV decay amplitudes are consistently under- 
estimated without annihilation contributions. When they are included, B + — > tt + K*° is 
dominated by fK'F^m^a^ + 6 3 (V, P) while B° — > K + p~ is dominated by f K AQ p m 2 B (a^ — 
r%a 6 )+b 3 (P,V) with 

b 3 (M 1 , M 2 ) = ^{CyLi(M 1; M 2 ) + C B Ai(M u M 2 ) + (C 5 + N e C 6 )Al{Mi, M 2 )} . (12) 



Since the penguin terms a 4 



—0.03 and a 4 — r^a% 

A, 



0.037 are of opposite sign, the 



key observation is that, the annihilation terms b$(y,P) and b%(P, V) are also roughly of 
the opposite sign, because numerically (C5 + N C C§)A{ is the dominant term in 63 while 
4(P, V) = -A{{V, P) if the annihilation parameter X A are the same for both channels. So 
with the inclusion of annihilation topology, QCDF can easily enhance both ratios without 
fine tuning. For example, using the default parameters given in ref. p] but letting the 
annihilation parameter = 1 we can see from Fig. 1 that the ratios are then consistent 
with experimental observations. We believe that this is a strong indication that annihilation 
topology probably plays an important role, at least in penguin-dominated B — > PV decays. 
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FIG. 1: The ratios B(B+ -> ir + K*°)/B(B° -» p+vr") (left plot) and B(B° -> K + p-)/B(B° -c 
/9~7r + ) (right plot) versus the weak annihilation phase <pA- The default parameters of ref. are 
used but letting the annihilation parameter pa = 1. The dashed lines show the ratios without 
weak annihilation contributions. The gray areas denote the experimental measurements with lex 
error. 



Similar disagreement can be found in the ratio of B° — > u K over B — > 7r + p~ . But we 
will not go into details here because B° — > uK° decay is quite similar to B° — > K + p~ decay. 
It may also be interesting to discuss B — > <pK decays, but again they are quite similar to 
B — > iiK* decays if there is no new physics contributions in this decay. 
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IV. DIRECT CP VIOLATIONS 



If one agrees that annihilation terms are very likely important in charmless B decays, one 
negative point is that, these intractable terms will bring large uncertainties to the theoretical 
predictions, not only in the branching fractions, but also in direct CP asymmetry . In the 
following we will try to see whether we can get some interesting relations with a few selective 
ratios of direct CP asymmetry (DCPVs). 

Let us first consider the ratio of DCPVs Acp(7r + 7r~) over Acp{it~ K + ). Belle has claimed 
large DCPV observed in B° — > 7r + 7r~ decay while BaBar has not confirmed it yet, but both 
of them are close to a measurement on Acp{^~ K + ) 19 



0.58 ±0.15 ±0.07 (Belle) , (-8.8 ± 3.5 ± 1.8)% (Belle) , 

A nn = < A n K = < (13) 

0.19 ± 0.19 ± 0.05 (BaBar) ; [ (-10.7 ± 4.1 ± 1.2)% (BaBar) . 

Since QCDF generally predicts small DCPVs for charmless B decays, at first glance it 
seems to contradict the Belle data badly. But notice that annihilation topology and other 
sources of power corrections could bring large uncertainties, it is of interest to have further 
investigations. In QCDF, we could define the tree and penguin amplitudes as: 

= UF B *m 2 B {a u x + a\ + r x a u 6 ) + f B fMh + h + 2b 4 ) = UF B *m 2 B T , 

= UF B *m 2 B {al + r x al) + fsUUh + 2b 4 ) = UF m m 2 B P , 

T nK = f K F B *m 2 B {a\ + a\ + r x <) + hfjxh = fKF B *m 2 B f , 

PnK = f K F B *m 2 B (al + r x a%) + f B Uf K h = f 7T F Bn m B P . (14) 

Then the DCPVs can be expressed as 

B (B° 7T+7T-) - B(B° -> 7T+7T-) A\V ub V ud V cb V cd T P\ sin 7 sin5 



A — 



B° -> 7T + 7T-) + B(B° -> 7T + 7T-) 2B(B -> TT + TT 



= B(B° -> -K+K-) - B(B° -> ti-K + ) A\V ub V us V cb V cs TP\sm 1 sm5 
* K ~ b ^ Ti+K-) + B(B° -> n-K+) 2B(B -> ir+K') ' 1 j 

where 5 = 5p — 5t is the strong phases difference between the penguin and tree amplitudes. 
It is easy to see that many factors cancel out for the ratio 



Anr 11 B(B -> ir+K- 



A„k f 2 K B(B - tt+tt- 



TP 



sin 5 . sin 5 . , 

— ~ (-2.7 ±0.3)-^ , (16) 
sin o sin o 



TP 

where the ratio TP/TP has been taken to be 1, which is a reasonable approximation in 
QCDF at about 10 percent level uncertainty. The experimental data on relevant branching 
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ratios |8( have been used in the above equation, and only the experimental uncertainties are 
included in the error estimation. 

At first sign, it is amazing to realize that one would expect very naturally a larger DCPV 
for 7c + tt~ decay compared with tc~K + decay, if the strong phases of both channels are not 
too small and of similar magnitude. Notice that the DCPV measurements already told us 
that at least the strong phase of 7r + 7r~ channel should not be small and QCDF predicts quite 
similar strong phases for both channels if other sources of power corrections are negligible. 
Actually, it is quite plausible that CP-violating asymmetry in 7r + 7r~ decay should be bigger 
than that for 7t~K + by a factor of 3 — 4 since the 7r + 7r~ decay rate is smaller than the 
7t~K + decay rate by a similar factor 3 — 4, given the fact that in QCDF the CP-violating 
parts of the tree-penguin interference terms are almost equal as shown above. The current 
experimental value |3| for this ratio is 

A n7T 0.42 ±0.13 , , 

-^L = = -4.2 ±1.8, 17 

A nK -0.10 ±0.03 ' v ' 

which is still consistent, within la error, with the theoretical estimation of —2.7 ±0.3 under 
the assumption that the strong phases 6 and 6 are the same. It should be interesting to 
keep an eye on this ratio and any significant deviation from theoretical estimation would 
suggest either different strong phases between tttt and irK decays or New physics effects 
affecting one of the channels. It is also interesting to note that inelastic B —>■ DD — > nir and 
B — > DD S — > Ktt rescatterings or charming penguin contributions could also produce a large 
strong phase, but since these contributions are related by SU(3) symmetry and CKM factor, 
the interference terms would be essentially equal and our relation for the CP asymmetry 



applies and we expect large CP asymmetry in B — ► as found recently 



21 



With the same reasoning, we can get similar relations for other decay channels such as: 



A CP (B° -> P+7T-) B(B° ^K*+7t-) f 2 p sin 5, 



A CP (B° -> K*+ir~) B(B° -> p+vr-) /£. sin 6 vK . ' 

A CP (B° - p-7T+) B(B° - p~K+) fl sin 6 p7T 



(18) 



Acp(B° ^ p-K+) B(B° ^ P-7T+) fl sin 8 pK ' 

If the above pairs of strong phases are roughly the same, which is true in QCDF, The DCPVs 
of the penguin-dominated decays would be about 1.5 times larger than those of their tree- 
dominated counter parts. Precise measurements on these ratios could help us to get some 
insight into the strong dynamics of B decays. 
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V. CONCLUSION 



In this paper, we have used the ratios of the branching fraction of B + — > ti + K* to that 
of B° — > 7r _ p + , and of 5° — * K + p~ to that of 5° — ► 7r~p + , to show clearly and with 
greatly reduced theoretical uncertainties, that power corrections in charmless B decays are 
probably large. The key observation is that QCDF predicts the annihilation terms for B + — > 
n + K*° and B° —>■ K + p~ are almost the same magnitude but opposite in sign. The result is 
that both these two branching fractions could be enhanced by the annihilation contribution 
to accommodate the experimental data without fine tuning. Assuming the presence of 
the annihilation contribution, we then derive relations between CP asymmetries for a few 
selective decay channels and shows that QCDF would predict the direct CP asymmetry of 
B 7r + 7r~ to be about 3 times larger than that of B — > tt with opposite sign, which 
is consistent, within one sigma error, to the current experimental data —4.2 ± 1.8. Any 
significant deviation from this prediction would suggest either new physics or possibly the 
importance of long-distance rescattering effects. We also discuss other similar direct CP 
ratios which might help us to get some insight into the strong dynamics of charmless B 
decays. 



Note Added 



After finishing this paper, we were informed that Eq. (|16j) . the relation between the 
direct CP asymmetries and CP-averaged branching ratios of B° — ► 7r + 7r~ and B° — > n~K + , 
has been derived previously by R. Fleischer 
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